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Abstract 

Spin correlations, using a generic spin basis, are investigated to leading order 
in QCD for top quark production at lepton colliders. Even though, these 
radiative corrections induce an anomalous 'y/Z magnetic moment for the top 
quarks and allow for single, real gluon emission, their effects on the top quark 
spin orientation are very small. The final results are that the top (or anti-top) 
quarks are produced in an essentially unique spin configuration in polarized 
lepton collisions even after including the 0{a s ) QCD corrections. 
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I. INTRODUCTION 



The discovery of the top quark, with a mass near 175 GeV |{]|§, provides us with a unique 
opportunity to better understand electro-weak symmetry breaking and to search for hints of 
physics beyond the standard model. It has been known for sometime that top quarks decay 
electroweakly before hadronization HQ and that there are significant angular correlations 
between the decay products of the top quark and the spin of the top quark p. Therefore 
if the production mechanism of the top quarks correlates the spins of the top and anti-top 
quarks, there will be a sizable angular correlations between all the particles, both incoming 
and outgoing, in these events. 

There are many papers on the angular correlations for top quark events |],[7j produced 
both at e + e~ colliders M and hadron colliders P,|TT| . In most of these works, the top quark 



spin is decomposed in the helicity basis. Recently, Mahlon and Parke |T0| have proposed a 
more optimal decomposition of the top quark spin which results in a large asymmetry at 
hadron colliders. Parke and Shadmi [TTJ extended this study to e + e~ annihilation process 
at the leading order in the perturbation theory and found that the "off- diagonal" basis is 
the most efficient decomposition of the top (anti-top) quark spin. In this spin basis the 
top quarks are produced in an essentially unique spin configuration. Since this result is of 
great interest, it is of crucial importance to estimate the radiative corrections to this process 
which are dominated by QCD effects. 

The QCD corrections to top quark production can be calculated perturbatively at ener- 
gies sufficiently above the production threshold of the top quark pairs. The analytical study 
of QCD radiative corrections to heavy quark production was pioneered in Ref. fl2| (see e.g. 



Ref. [13] for a recent article). Polarized heavy quark production, in the helicity basis, has 



also been investigated by many authors 

In this article, we present an analytic differential cross section for polarized top quark 
production at the QCD one-loop level. We focus on the issue of what is the optimal de- 
composition of the top quark spin for e + e~ colliders Q We have calculated the cross section 
in a "generic" spin basis which includes the helicity basis as a special case. The radiative 
corrections, in general, add two effects to the tree level analysis: the first is that a new 
vertex structure (anomalous 7/Z magnetic moment) is induced by the loop corrections to 
the tree level vertex, the second is that a (hard) real gluon emission from the final quarks 
can flip the spin and change the momentum of the parent quarks. Therefore, compared to 
the radiative corrections to physical quantities which are spin independent, it is possible 
that spin-dependent quantities maybe particularly sensitive to the effects of QCD radiative 
corrections. 

The article is organized as follows. In Section 2, we examine the QCD corrections to 
the polarized top (anti-top) quark production in the soft gluon approximation. The aim of 
this section is: (1) we estimate the numerical effects from the new vertex structure on the 
spin correlation found in the tree level analysis and (2) we show that we can use the off- 
diagonal basis as a optimal basis also at the QCD one- loop level. In Section 3, we present our 



1 The physics of top quark production at muon colliders and e + e colliders is identical provided 
the energy is not tuned to the Higgs boson resonance. 
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analytic calculations of the full one-loop corrections to the polarized top quark production 
in a generic spin basis. We give the numerical results both in the helicity, beamline and 
the off-diagonal bases in Section 4. Here we compare the full one loop results with those 
of the tree level and soft gluon approximations. Finally, Section 5 contains the conclusions. 
The phase space integrals which are needed in Section 4 are summarized in Appendix A. 
The unpolarized total cross section for top pair production, using our results, is given in 
Appendix B as a cross check. 



II. SPIN CORRELATIONS IN THE SOFT GLUON APPROXIMATION 



In this section we derive the first order QCD corrected spin dependent, differential cross 
section for top quark pair production in the soft gluon approximation(SGA). It is instructive 
to first consider the soft gluon approximation because in this approximation only the QCD 
vertex corrections modify the spin correlations of the top quarks. The full one-loop analysis 
will be given in the next section. We use the same generic spin basis as in Ref. ||11|| . In this 
paper we do not consider transverse polarization of the top quarks f\ since we are interested 
in how QCD corrections modify the tree level spin correlations and which spin basis is the 
most effective for spin correlation studies. Therefore we use a generic spin basis with the 
spin of the top quark and anti-top quark in the production plane. We define the spins of 
the top and anti-top quarks by the parameter £ as given in Fig.l. The top quark spin is 
decomposed along the direction s t in the rest frame of the top quark which makes an angle £ 
with the anti-top quark momentum in the clockwise direction. Similarly, the anti-top quark 
spin states are defined in the anti-top rest frame along the direction Sf having the same 
angle £ from the direction of the top quark momentum. We use the following notation in 
this paper: the state if (ti ijj refers to a top with spin in the +s t (— s t ) direction in the 
top rest frame and an anti-top with spin +sj (— sj) in the anti-top rest frame. 

The one-loop QCD correction to the cross section is given by the interference between 
the tree and one-loop vertex diagrams in Fig. 2. At the one- loop level, the 7 — t — i and 
Z — t — i vertex functions can be written in terms of three form factors A, B, C as follows: 
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where Q t = f is the electric charge of the top quark in units of the electron charge e, 9w is 



2 It is known that the transverse top quark polarization becomes nonzero when the higher order 
QCD corrections are included and that this transverse polarization is very important and related 
to the phenomena of CP violation. 
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the Weinberg angle, and m and (t M ) are the mass and the momentum of the top (anti-top) 
quark (7^/ L = 7 M ^^). The top quark couplings to the Z boson are given by 



Qt — 



3-4 sin 2 9 W 



2 sin 2 9 W 



Q R = 

6 cos 9w 3 cos 6 W 



(3) 



After multiplying the wave function renormalization factor (we employ the on-shell renor- 
malization scheme), the "renormalized" form factors read 
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where (3 is the speed of the produced top (anti-top) quark and the strong coupling constant 



is a s 



c 2 (R) 

4tt 



= u^J g 2 with C2(R) = 3 for SU(3) of color. We have introduced an 
infinitesimal mass A for the gluon to avoid infrared singularities. In the above expressions, 
we have shown only the real part of the form factors because (1) the Z width is negligible 
in the region of center-of-mass (CM) energy >/i far above the production threshold for top 
quarks and (2) we are not considering the transverse polarization for the top quarks. The 
contribution from C can be neglected since it is proportional to the electron mass. 
The differential cross section at the one loop level is given by 
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x (A Lfl sinf + fl LJi cosf ±£> Lfl )(l + 2A + 2B) 
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Here, the angle is the scattering angle of the top quark with respect to the electron in 
the zero momentum frame, a is the QED fine structure constant and 7 = 1/ y/1 — (3 2 . The 
quantities A LR , B LR , B LR , D LR and D LR are defined by 
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with 



//j = -Q t + QlQ 



l 



e "** sin 2 0vK s-Mf ' 

where Mz is the Z mass (as mentioned before, we neglect the Z width) and I, J G (L, R). 
The electron couplings to the Z boson are 

r _ 2 sin 2 9 W - 1 B _ sin 2 6 W 

2 COS t>w cos o'vk 

The cross sections Eqs.(|7],^) contain an infrared singularity (in the form factor A) that 
will be canceled by the contributions from the real gluon emission. In the soft gluon ap- 
proximation, it is very easy to calculate the real gluon contribution. As is well known, the 
amplitude for the soft gluon emissions can be written in the factorized form proportional 
to the tree amplitude. This means that the soft gluon emission does not change the spin 
configurations or momentum of the produced heavy quark pairs from the tree level val- 
ues. Therefore the QCD radiative corrections enter mainly through the modifications of the 
vertex parts Eqs.([l],Q). The cross section for the soft gluon emissions can be written as 

_da__ da 
a cos v a cos U 

where the subscript denotes the tree level cross section. The soft gluon contribution Jjr 
is defined by 

/■fc°=a; max rl^h ft f \ 2 

J m = -4irC 2 (R)a s 



where w max 
we obtain 
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(2ir) 3 2k \t ■ k t-kj ' 
is the cut-off of the soft gluon energy. This integral can be easily performed and 
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By adding the one loop contributions Eqs.(0J|) and the soft gluon ones Eq. fllOD , one can 
see that the infrared singularities, In A, are canceled out and the finite results are obtained 
by replacing 2A by 
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in Eqs.(0,|). 

The cross sections for e^e\ can be obtained by interchanging L, R as well as f, j in the 
above formulae. 

Since we are interested in maximizing the spin correlations of the top quark pairs we 
must vary the spin angle, £, to find the appropriate spin basis. At tree level, it is known 
that there exists the "off-diagonal" basis which makes the contributions from the like-spin 
configuration vanish [IT|. At order 0(a s ), we find that definition of the off-diagonal basis 
for e^e^ scattering is not modified by QCD corrections, with the spin angle, £, satisfying 
the tree level relationship 
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The first order QCD corrected cross sections in this basis are 
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A similar result holds for e^e\ scattering. 

In Fig. 3 we show the differential cross sections in the off-diagonal basis, Eq. (|TTD , for 
yfs = 400 GeV. The following values for the parameters of the standard model were used 



m 



My 



top 



175 GeV 



M, 



91.187 GeV 



a = — , a s (M 2 z ) = 0.118 , sin 2 9 W = 0.2315 . 
128 

We use y/s as the renormalization scale and the pole mass for the top quark. For e^e^ 
scattering the up- up (tfff) and the down-down tfii components are identically zero. The 
total cross section is more than 99% up-down tft^ and less than 1% down-up tfa. For e^e\ 
scattering the up-up (t-pt^ ) and the down-down t^i components are non-zero because we 
have used the off-diagonal basis for e^e^ scattering. However the down- up t\t^ component 
is still more than 99% of the total cross section. 

Although there exists a magnetic moment modification to the 7 / Z — t — i vertex from 
QCD corrections, this does not change the behavior of the spin dependent cross sections in 
the off-diagonal basis. The QCD corrections, however, make the differential cross sections 
larger by ~ 30% compared to the tree level ones at this y/s. Thus the off-diagonal basis 
continues to display very strong spin correlations for the top quark pairs even after taking 
the QCD corrections into account, at least in the soft gluon approximation. 
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In the previous paragraph the cut-off energy for the soft gluon has been chosen as uj max = 
10 GeV. The results, of course, depend on the value of u; max . The u; max dependence of the 
cross section is examined in Fig.4. The cross section behaves quite uniformly as the value of 
w max is changed thus the above conclusions remain qualitatively the same for any reasonable 
value of u; max . 



III. SINGLE SPIN CORRELATIONS IN E+E~ PROCESS 

The soft gluon approximation used in the last section has two shortcomings. First, soft 
gluon emission cannot change the spin of the heavy quarks and second the heavy quark 
pairs are always back to back. Neither of these approximations is valid for hard gluon 
emission. Hence it is possible that the full 0(a s ) QCD corrections might completely change 
the conclusions of the previous section. Therefore, in this section, we investigate the full 
0(a s ) QCD corrections. Since, in the presence of a hard gluon, the top and anti-top quarks 
are not generally produced back-to-back, it is more sensible to consider the single heavy 
quark spin correlations, namely the inclusive cross section for the production of the top (or 
anti-top) quark in a particular spin configuration. We have organized this section as follows. 
After defining the kinematics and our conventions, we give the polarized cross section for 
the top quark using a generic spin basis closely related to the spin basis of the previous 
section. The numerical analysis will be relegated to the next section. 



A. Amplitudes and Kinematics 



The principle result of this section will be the inclusive cross section for polarized top 
quark production, e + e~ — > or t[ + X, where X = i or ig. (The cross section for the 
anti-top quark inclusive production can be easily obtained from the results in this section.) 
Since the vertex corrections are the same as in the previous section all that is required is the 
full real gluon emission contributions. The real gluon emission diagrams to leading order 
in a s are given in Fig. 5. This figure also defines the momenta of particles. We will use the 
spinor helicity method for massive fermions |10[ to calculate the squares of these amplitudes 
for a polarized top quark. The top quark momentum t is decomposed into a sum of two 
massless momenta t\ , t 2 such that in the rest frame of the top quark the spatial momentum 
of ti defines the spin axis for the top quark. 
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where s t is the spin four vector of the top quark. 
The amplitude for Fig. 5 is given by 
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where a,f3 are the spin indices for quarks and e is the polarization vector of the gluon. T a 
is the color matrix. The coupling constants an are defined as follows: 

o-li _ e 2 g 

The expressions for the squares of the amplitudes given below have been summed over 
the spins of the unobserved particles (the anti-top quark and gluon) as well as the colors of 
the final state particles. Let us write the square of the amplitude for the top quark with 
spin "up" as 
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where u± = and all momentum, p, under the "Tr" operator are understood to be ^. 
By interchanging the t x and t 2 vectors in the above expressions, we can get the amplitude 
square for the top quark with spin "down". Since we neglect the Z width in this paper, all 
the coupling constants an are real. 

To define the spin basis for the top quark we naturally extend the spin definition of the 
previous section to the present case. The top quark spin is decomposed along the direction 
s t in the rest frame of the top quark which makes an angle £ with the sum of the anti-top 
quark and the gluon momenta in the clockwise direction, see Fig. 6. 

To calculate the cross section from Eq.fli~5D, we take the CM frame in which the e + e~ 
beam line coincides with the z-axis, 
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9=^(1, 0,0,1) , 9=^(1,0,0,-1). 

We specify the variables x, y and z which are related to CM energies of the gluon, top and 
anti-top quarks by 

= 1 2k-(q + q) = 2t-(q + q) = 2i ■ (q + q) 

x — i , y — i , z — i 

s s s 

The momenta of the final state particles, in terms of these variables, are 



k = ^(l-x, (l-x)k) ,t = y^(l-y,a(y)i) , t = ^(1 - z , a(z) t) , 
where A means the unit space vector and 



a(y) = yj(l- y y-a , a(z) = yftl 



z) 2 



with a = Am 2 /s. Fig. 7 defines the orientation of the top and anti-top momenta and by 
energy-momentum conservation the momentum of the gluon is also determined.^ One can 
easily obtain the spin four vector of the top quark in the CM frame by boosting the spin 
vector characterized by £ in the top quark rest frame in the direction of top quark momentum 
by (3(y) (the speed of the top quark in the CM frame). The explicit form for t\ (t 2 = t — t\) 
is given by 

*? = j[7(v)(l-0(v)cosfl] , 

t{ = — [7(2/) (/3 (y) — cos£) sin^cos^ + sin £ cos # cos ] , (16) 

t\ = — [j(y) ((3(y) — cos£) sin 6 sin ip + sin£ cos 9 sin ip] , 
o m 



where 



h = y [ (/%) - cos cos 6 1 - sin £ sin ] 



aj(y) = l-y , y/aj(y) (3(y) = a(y) . 



If we eliminate the gluon momentum k using the energy momentum conservation and 
use the angular variables x > i n Fig- 7 to specify the orientation of the anti-top quark (if 
one eliminates the anti-top momentum, one can proceed in the similar way by introducing 
other angular variables), the square of the amplitude Eq . (|15|) can be written as 



\M(e~ L e + R -> t^tg)\ 2 = J N C C 2 (R) a\ L M x + a 2 LR M 2 + a LL a LR M s , (17) 



3 Our Fig. 7 contains an extra degree of freedom, 
These works of Tung et al. correspond to setting 
unable to reproduce their results. Our variable (ft 
of ref. [@. 



(ft, compared to Fig. 3 of the first paper in ref. JT^] . 
our (ft = 7r/2. Due to this difference we have been 
corresponds to the variable <j>i2 in the first paper 
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where M, are the functions of y , z, angles defined above and the spin orientation £. 
2 r. 
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+ — 

yz 
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In the above equations, 9k is the angle between the z axis and the gluon momentum 

(y + z) cos 9k = —a(y) cos 9 — a(z) cos 9 , 

and St is defined as 

St EE -(*! -t 2 ) . 

S 

Using Eq. (|l^) we find that the products of St with momenta q , q and t are 

St ■ q = {(1 — 2/) cos 6* — a(y)} cos£ + v^sin^sin^ , 
St - q — { — (1 — y) cos9 — a(y)}cos^ — ^/asm.9 sm.^ , 
St -t = { — (1 — z)a(y) + (1 — y)a(z) cosx} cos£ 

+ y/aa(z) sin £ sin x cos (j) . 
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The unpolarized top quark production process is given by dropping the spin dependent parts 



(terms proportional to 5t) in Eq. (|17|) . As a check we have reproduced the results of Ref. [12 
by putting a LL = a LR = ~^Q q - 
The cross section is given by 
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where (PS) 3 is the three particle phase space. 
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We introduce a small mass A for the gluon to regularize the infrared singularities. It is easy 
to rewrite the above phase space integral as, 
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where dQ = d cos ddip is the solid angle for the top quark and a\ = X 2 /s. The integration 
regions over y and z are determined by the condition | cosxl < 1? 

y + = 1 - y/a , y_ = ^/aa^ + a x , 

= Ay + a y ( X ~ ^ ~ I + aA ) + aA ± a ( y ^( y ~ ax ^ ~ aax 
The integration over the angle is not difficult if one uses the relation 

cos 9 = cos 6 cos x + sin sin x cos . 
The integrals we need are the following: 
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Due to the 5 function in the phase space integral, the angle x is a function of y and z: 
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where we have put a\ to be zero since the A — > limit does not produce any singularities 
in the (squared) amplitude. The remaining integrals to get the cross section are over the 
variables y and z. According to the type of integrand, we group the phase space integrals 
(after the integrations over the angular variables) into four distinct classes { J{\ , {A^} , {Li} 
and {Ki} ||15|| . The individual integrals of these classes are summarized in Appendix A. 



B. Cross Section in the Generic Spin Basis 



We write the inclusive cross section for the top quark in the following form. 

( e L 4 h X ) = ^TT H {Dklmn + OisCklmn) COS* 9 SUl' 9 COS™ £ SUl™ £ 



dcos9 



As 



(19) 



klmn 



where D kimn are the contributions from the tree and the one-loop diag rams and Ckimn &re 
from the real emission diagrams. Let us first write down the D/,i mn , 



-Doooo 

-D2000 
-D1000 
-D0010 
-D2010 
-D1010 

-D0101 
-D1101 



/3[fl L + fl R + 2af LL f LR ](l + a s V I ) 

-PWll + Ilr) 2 - PVll - fL R f]&sVn 

+ + &sVj) + PHfLL - fLRf&sVu , 

2(3 2 (f 2 LL -fl R )(l + a s V I ) , 

Wll + Ilr) 2 + [3 2 (f L L - /lk) 2 ](1 + & s Vj) 

-2(3[(f LL + f LR f - p 2 (f LL - f LR ) 2 ]a s V n 

-^{/ll + f L R) 2 [a(l + a s V!) - (1 + a)a s V n ] , 



^(f 2 LL - f 2 LR )[a(l + a s Vj) - (1 - a)a s V n } , 



with 



(3 = (3(0) = VT 

A , B are defined in Eqs. djj^) . 
For the C^ mn we find, 



a s Vi = 2A + 2B 



a.V n = B . 



a 



0000 



a 



2000 



2 (/El + fin) 
+ 4a f LL f LR 

2 (fL + fl R ) 
+ 4a f LL f LR 



Jjk - (4 - a) J 3 + (2 + a) J 2 + a.J x + 



^IR — 4 J 3 — J 2 — J\ + ^-^2 



;i - a)Jl K - (4 - a) J 3 + (2 - a) J 2 - aJ x - 



<^2 + J\ — 



C1000 = 2 (/| L - /E B ) (1 - a)J 2 K + aiV 10 - 2(4 - 3a)N 9 
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(4 - 5a)iV 8 - 2N 7 + 2N G + 6N 3 + 2N 2 



Cooio 

C2010 
C1010 



C 



1000 



ill fill) 



AN 6 - aN A - aN 3 



1 



aN 2 + (4 - a)iVi + -R 3 



+ (fl 

Cooio - 2 (/ 2 L ~ /Ei?)-^3 ) 

UIl + flu) [ 2(2 - a) A + 2aJ 4 - 2(8 - 5a) J 3 

- 2a(l - a)L 8 - 2a(l - a)L 7 - 2(4 + 3a)L 6 + 2(4 - 3a)L 5 
+ 2aL 4 + a(10 - a)L 3 + (8 - 6a - a 2 )L 2 - 2(4 - 3a + a 2 )L 4 



+ 2a f LL f LR 



2^ir 



CI 



0101 



(/ll + fill) 



Mr 



^ J 3 + 4L 6 + 4L 5 + aL 3 + aL 2 - 2(2 - a)LiJ , 

(8 + a) J 3 - (8 - a)(l - a)L 7 - (12 + a)L 6 
- 2aL 5 + 2aL 4 + (8 + a)L 3 + (4 - 5a)L 2 + 4(2 - a)L x 
+ fLLfLR [ 4Jm - (16 - a) J 3 - a(l - a)L 7 + (4 + a)L 6 
+ 2aL 5 + aL 3 - (4 - 3a)L 2 + 2aL 1 ] , 



C 



1101 



C1000 + \/a (/ll - 



X 



i(12 + a)iV 3 + ^iV 2 -(12 + a)iV 1 



-aiV 10 - 9aiV 9 + 2iV 7 - 27V 6 

- a(l + a) #8 + a(l - cl)K 7 + 9a(l - a)ir 6 - 3a(5 + a)K s 

- (12 + a)K A + (4 - 5a)K 3 - 3(4 + ha)K 2 + 3(4 + a)(l - a)K x 



where we have defined 



J IR 
J 2 

Ri 



(2 - a)J 6 - aJ 5 , 
2(2 - a)N 13 - aN 12 - aN u , 
-4a(l - a)L 7 - 4(2 - a)L 6 - 8(1 - a)L 5 
+ a 2 L 4 + a(2 + 3a)L 3 - a(2 - a)L 2 + (8 - 8a + 3a 2 )Li , 
R 2 = -4L 6 - 4L 5 - aL 3 - aL 2 + 2(2 - d)L\ , 
R 3 = a 2 N 10 + 3a(2 + a)N 9 + AaN 3 - 2aN 2 + 8(1 + a)iVi 

+ 2a 2 if 8 - a 2 (l - a)K 7 - 3a(l - a) (2 + a)K 6 + 6a(l + 2a)K 5 
+ 2(4 + 3a)fT 4 + (a 2 + 6a - 8)if 3 + 3a(8 + a)K 2 - 12a(l - a)K x . 

Note that the integrals J} K , J 2 R namely J5 , J6 , AT n , N X2 and iVi 3 contain the infrared sin- 
gularity. This singularity is exactly canceled out in the sum Eq. (|19|) by the contributions 
from D k i mn . We are now ready to discuss our numerical results. 
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IV. NUMERICAL RESULTS 



We are now in the position to give the cross section for polarized top quark production 
for any e + e~ collider. Since we did not specify the spin angle £ for the top quark, we can 
predict the polarized cross section for any top quark spin. The spin configurations we will 
display are the helicity, the beamline and the off-diagonal bases for center of mass energies 

= 400 GeV , 800 GeV , 1500 GeV. 

Table 1 contains the values of the maximum center of mass speed, the running a s , the 
tree level cross section, the next to leading order cross section and the fractional 0(a s ) 
enhancement of the tree level cross section (k) for top quark pair production in e^e + scat- 
tering. 





400 GeV 


800 GeV 


1500 GeV 


p 


0.484 


0.899 


0.972 


a s {s) 


0.0980 


0.0910 


0.0854 


OTotai Tree (pb) 


0.8708 


0.3531 


0.1047 


a Total 0{a s ) (pb) 


1.113 


0.3734 


0.1084 


K 


0.2783 


0.05749 


0.03534 



Table 1: The values of (3, a s , tree level and next to leading order cross sections and k 
for e^e + scattering. 

At s/s = 400 GeV the QCD corrections enhance the total cross section by ~ 30% compared 
to the tree level results whereas at higher energies, 800 and 1500 GeV, the enhancements 
are at the ~ 5% level. 

First we will show the numerical values for the coefficients {Dkimn + ^sCkimn) in Eq. (|l^) . 
Since the dominant effect of the 0{a s ) corrections is a multiplicative enhancement of the 
tree level result we have chosen to write 

Dklmn + C^sCklmn = (1 + K )-Dfc£ mn + Sklmn- 

The (1 + K)D® lmn terms are the multiplicative enhancement of the tree level result whereas 
the Skimn give the 0(a s ) deviations to the spin correlations. The numerical value of these 
coefficients are given in Table 2. 
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/ — 

Vs 


400 GeV 


800 GeV 


1 V f~\f~\ f~ \ T T 

1500 GeV 


(1 + k)D» 000 


1.511 


1.722 


1.671 


(1 + k)D% 00 


0.2404 


1.241 


1.526 


(1 + k)D% 00 


0.7809 


2.126 


2.406 


(1 + k)D° 00W 


0.3905 


1.063 


1.203 


(1 + k)D% 010 


0.3905 


1.063 


1.203 


(1 + k)D% 10 


1.751 


2.963 


3.197 


(1 + k)D° 101 


1.452 


1.100 


0.6186 


(1 + «)DJ 101 


0.3416 


0.4650 


0.2807 


5*0000 


-0.002552 


0.0005645 


0.0H72 


5*2000 


0.007655 


-0.001682 


-0.03516 


5iooo 


0.02154 


0.01224 


-0.02085 


5ooio 


0.01094 


0.01586 


0.008357 


^2010 


0.01059 


-0.006158 


-0.03614 


5*1010 


0.004433 


-0.02030 


-0.06134 


5oioi 


-0.006564 


-0.04413 


-0.04952 


5noi 


0.007920 


-0.01270 


-0.02047 



Table 2: The values of (1 + K)D% lmn and Ski mn for e L e + scattering. 



The ratios Skimn/ 0- + K )Dli mn are never larger than 10% and are typically of order a few 
percent. Hence the 0(a s ) corrections make only small changes to the spin orientation of the 
top quark. 

To illustrate the different spin bases we present the top quark production cross section 
in the three different spin bases discussed in ref. [ll|. One is the usual helicity basis which 
corresponds to cos£ = +1. The second is the beamline basis in which the top quark spin is 
aligned with the positron in the top quark rest frame (£ for the beamline basis is obtained 
from Eq.flllD with Jlr = 0). The third corresponds to the off-diagonal basis which has 
been defined in Eq.([ll]). Note that as /? — ► 1, all of these bases coincide. Therefore, at an 
extremely high energy collider, there will be no significant difference between these bases. 

In Fig. 8 we give the results for y/s = 400 GeV for both e~[e + and e^e + scattering using 
the helicity, beamline and off-diagonal spin bases. This figure shows the tree level, SGA and 
the full QCD results for all three spin bases. Since the results in the SGA almost coincide 
with the full QCD results the probability that hard gluon emission flips the spin of the top 
quark is very small. Clearly, the qualitative features of the cross sections remain the same 
as those in the leading order analysis. That is the top quarks are produced with very high 
polarization in polarized e + e~ scattering. In Table 3 we give the fraction of the top quarks 
in the sub-dominant spin configuration for ej / e + scattering, 



e T e 



t^iiorig) K 



Total 



for the three bases. Similar results also hold for e R e 



scattering. 
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v/i = 400 GeV 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.336 


0.0119 


0.00124 


SGA 


0.332 


0.0113 


0.00129 


0(a s ) 


0.332 


0.0115 


0.00150 



Table 3: The fraction of the e^e + cross section in the sub-dominant spin at 
■y/s = 400 GeV for the helicity, beamline and off-diagonal bases. 
For the soft gluon approximation (SGA) u max = (y/s — 2m)/ 5 = 10 GeV f\ 

In Fig. 9 and 10 we have plotted the similar results for 800 and 1500 GeV colliders. The 
fraction of top quarks in the sub-dominant spin component for ej j e + is given in Tables 4 and 
5. 



y/s = 800 GeV 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.168 


0.0690 


0.0265 


SGA 


0.164 


0.0679 


0.0272 


0(a s ) 


0.165 


0.0708 


0.0319 



Table 4: The fraction of the ej j e + cross section in the sub-dominant spin at 
y/s = 800 GeV for the helicity, beamline and off-diagonal bases. 
For the soft gluon approximation (SGA) u> max = (y/s — 2m)/ 5 = 90 GeV. 



= 1500 GeV 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.132 


0.0978 


0.0466 


SGA 


0.130 


0.0973 


0.0472 


0(a s ) 


0.133 


0.101 


0.0552 



Table 5: The fraction of the ej j e + cross section in the sub-dominant spin at 
y/s = 1500 GeV for the helicity, beamline and off-diagonal bases. 
For the soft gluon approximation (SGA) u max = (y/s — 2m)/ 5 = 230 GeV. 

Our numerical studies demonstrate that the QCD corrections have a small effect on 
the spin configuration of the produced top (or anti-top) quark for any spin basis. The 
off-diagonal and beamline bases are clearly more sensitive to the radiative corrections than 
the helicity basis. However, in the off-diagonal basis, the top (and/or anti-top) quarks are 
produced in an essentially unique spin configuration even after including the lowest order 
QCD corrections. 

V. CONCLUSION 

We have studied the 0(a s ) QCD corrections to top quark production in a generic spin 
basis. The QCD corrections introduce two effects not included in the tree level approxima- 
tion. One is the modification of the coupling of the top and anti-top quarks to 7 and Z 



4 In the SGA the fractions in Tables 3, 4 and 5 have a very small dependence on uj- 
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bosons due to the virtual corrections. The another is the presence of the real gluon emission 
process. First, we consider the QCD corrections in the soft gluon approximation to see the 
effects of the modified 7/Z — t — t vertex. Using this approximation, we find the tree level 
off-diagonal basis continues to make the like spin components vanish and that the effects of 
the included anomalous magnetic moment are small. 

Next we analyzed the full QCD corrections at one loop level. When we consider the three 
particle final state, the top and anti-top quarks are not necessarily produced back to back. 
So we have calculated the inclusive top (anti-top) quark production. In this paper we have 
given an exact analytic form for the differential cross section with an arbitrary orientation 
of the top quark spin. 

Our numerical studies show that the 0(a s ) QCD corrections enhance the tree level result 
and only slightly modifies the spin orientation of the produced top quark. In the kinematical 
region where the emitted gluon has small energy, it is natural to expect that the real gluon 
emission effects introduce only a multiplicative correction to the tree level result. Therefore 
only "hard" gluon emission could possibly modify the top quark spin orientation. What 
we have found, by explicit calculation, is that this effect is numerically very small. The 
size of the QCD corrections to the total cross section and the enhancement of the tree level 
results can be read off from the values of k in Table 1 of Sec. 4. At ^/s = 400 GeV, the 
enhancement is ~ 30% whereas at higher energies, 800 and 1500 GeV, it is at the ~ 5% level. 
Near the threshold, the QCD corrections have a singular behavior in (3, the speed of the 
produced quark, this factor enhances the value of the correction at smaller energy. The size 
of these corrections is reasonable for QCD, on the other hand, the change of the orientation 
of top quark spin are quite small. The deviation from the enhanced tree level result is less 
than a few percent. We can, therefore, conclude that the results of the tree level analysis 
are not changed even after including QCD radiative corrections except for a multiplicative 
enhancement. This means that for the beamline and off-diagonal bases, the top (and/or 
anti-top) quarks are produced in an essentially unique spin configuration. Actually, the 
fraction of the top quarks in the dominant (up) spin configuration for e^e scattering is more 
than 94% at all energies we have considered. 

As has been discussed in many articles, there are strong correlations between the ori- 
entation of the spin of the produced top (anti-top) quark and the angular distribution of 
its decay products. Therefore, measuring the top quark spin orientation will give us impor- 
tant information on the top quark sector of the Standard Model as well as possible physics 
beyond the Standard Model. 
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VII. APPENDIX A : PHASE SPACE INTEGRALS OVER Y AND Z 



The phase space integrals necessary to derive the cross section are summarized in this 
Appendix. Although many of them have already appeared in the literatures |I3j,15|, we 
will list all of them below for the convenience of the reader. After the integration over the 
angular variables, we are left with the following four types of integrals: 

Ji = f dydzfi(y,z) , Ni = [ , == =fi(y,z) , 

dydz f dydz 

-fi{y,z) , Ki= j- ^ T^fi{y,z) ■ 



{l-yY-a J ^ " 1 J {(1 - „)2 _ a }s/2 



The infrared divergences are regularized by the small gluon mass A and 8 = yl 



1 — Am 2 / s. For the type L{ integrals, a shorthand notation oo = J (1 — \fa)/ (1 + y/a) is 
used. The Ki integrals have a spurious singularity at the upper bound of the y integral, 
y + = 1 — yfa. Since this singularity turns out to be canceled out in the cross section, we 
regularize each integrals by deforming the integration region as y + —>■ 1 — \fa — e . Li 2 
is the Spence function. 
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N- 



10 



N- 



11 



N- 



12 



N- 



13 



dz 




y_ _ 4 


Va-^) 2 - 


a 


z 2 a 






1 


v^-^) 2 - 


a 


y 2 


r a 2 

— In h In a 

s 


+ 21n(l 


dy dz 




1 


^~v) 2 - 


a 


I 2 " 


A 2 

— In In a 

s 


+ 21n(l 


dy dz 




1 


xj(i-y) 2 - 


a 






In 


* + i 

s 2 







In 



+ 21n2 



1-/^ 
.1 + 0, 



1 



• ,2 fUL) .I / yg-Vg \ / 2^(1 - V5) ^ 
+ 8 \y/H(2-y/a)) + 8 {"'2 ) {(1 - y/E ~ P) 2 , 



2 

+ p 

1 

+ P 



Ll2 1 (1+/3) 2 r Ll2 



.1 + 0, 



-Ll2 « + Ll2 



Class -L Integrals: 



^2 



/ 



dy dz 1 — a , fl + 8 s 
2 In 



= -1 + 



(l-y) 2 -a 

z 

(1 - y) 2 - a y 

12 24 



4-a Xl-P, 



4-a (4-a) 2 , 
dyd^ y 



In 



'1 + 8' 



2/3 
4-a 



(1 — y) 2 — a z 



2 U-0, 



- In a + ln(2 + y/a) -\n(l + y/a) — 2 In 2 



+ 



+ 



I - y/a 



Li 2 (a;) + Li 2 (^5" 



-a;) 



Li 2 (i±^) + Li 2 ((2 + ^)i±^) + Li 



1 + w s 



2v^ 



i (2 + Va)(l + w) i 
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u 

L 6 



dy dz y 2 



(1 


-y) 2 - 


a 




dy dz 




(1 


-y) 2 - 


a 




dy dz 




(1 


-y) 2 - 


a 




dy dz 




(1 


-y) 2 - 


a 




dy dz 




(1 


-y) 2 - 


a 



In 



.1-0, 



2(3 



1 

l + -a 



3a 



6 



4 - a 
2 



4-a 
1 

la 

a \l-(3 



4 - a 
Li 2 (cu) + Li 2 



In 
In 



2 + ^/a 



+ 

.1-0, 

l : 0, 
a ^ 



+ 



a 



2(3 
4-a 

(uj — > — a;) 



Class _K~ Integrals: 



# 2 



K A 



dy dz 



{( l_ y) 2_ a} 3/2 



2(1 -y/a) 2y/E{l-y/a) 

y 



V^(2 - Va~) 2 
dydz 



4a ln (l_+ y/a){2- ^E) 2 



(4-a) 



{( l_ y)2 _ a} 3/2 



2(1 -v^) 2 h Va(Wa) + 



v^(2 - ^) 2 
2v^ 



e(l + y/a) (4-a) 2 



In 



2y / aa 



(l + v^)(2- y^) ; 
2y / aa) 



AT 3 = 



+ In 



(1 + V^) 

dz 



{( l_ y) 2_ a} 3/2 

2(l-v^) 2 2^/0(1-^/0) a(a 2 + 20a -32) (l + v^)(2-^) ; 
(2 - v^)3 11 e(l + v^) + 2(4 -a)3 11 

1 2^/a 2a(l - ^/a) 

~2 n (l + v^) + (4-a)(2- ^) 2 
dy dz 



2y / aa 



{( l_ y) 2_ a} 3/2 



2(1 - v^) 3 ^ Vgl - y/q) a 2 (12-7a) ^ (1 + ^)(2 - v^) 2 



(2 - v^) 3 e(l + Va) 2(4 - a) 



2^/aa 



1 , 2^ 2 

■In ,. v — + 



2 (1 + v/a) (4 - a) 
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{(l-y) 2 -a} 3 / 2 I 
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2(1 - y/a) ln Vg - y/g) (4 - 3a) ^ (1 + ^(2-^) 
a(2 - ^/a) e(l + ^/a) a(4 - a) 

4(1 - v^) 

a 



2y / aa 



a (1 + v/a) a(2 - ^a) 



+ 



dy dz 



1 

{{I - y) 2 - a}^ I 
2 



ln 2y^(l - y/a) (4 + a) ^ y/a)(2 - y/a) 



a(2 - y/a) e(l + ^a) 
2^ 4 



a(4 — a) 



2 v/a. 



aa 



-it 



+ 



/ 



a (1 + v/a) a(2 - y/a) a/3 
dy dz y 



2 2y/E(l-y/0) 



{(l_ y )2_ a }3/2 Z 2 a ^ € (l + y/a) 
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{(1 - y) 2 - a} 3 / 2 2 2 
2(1 - y^) ln /V^Wa)' 
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VIII. APPENDIX B : TOTAL CROSS SECTION 



Let us examine the unpolarized total cross section including the top quark pairs using our 
formulae. It is given by 

a T (e'e + ->t + t + X) = - [(y T {e~ L e^ -> t + i + X) + o T {e R e + L ^t + i+X) 



Note that only k = 0, 2 and I, m, n = terms in Eq . (|T9|) contribute to the total cross section. 
Integrating over the angle 8, we get 



°r(e£4 -> t + t + x ) 

>,2 



ira 



Ull + /l«) 2 (3 - P 2 ) 
6 



x 11 + & s jVr - 
+ 2(f LL -f LR ) 2 p 2 



V„ + -J 1 



Jo 



X (l + 6, ^Vj + 2Vfr + | J] R - | J a + 2(3 /?3 /?2) ^2 + ^ Jij 

The cross section for the process e^e^ is obtained by interchanging the coupling constant 
L <-> i? in the above expression. Parameterizing the total cross section as 



i? t ( S ) = ^J-j a T (e-e + 
>(0)/ „\ , 



7T 



tf 2 (J2)J# ) ( fl ) + 



where <7 pt = 47ra 2 /3s, we get the following numerical results at the CM energies ^Js 



400 , 500 , 1500 GeV. 



(GeV) 


R?\s) 


C 2 (R)R { t 1) (s) 


400 


1.0083 


8.9963 


500 


1.4190 


6.0267 


1500 


1.7714 


2.2911 



These are consistent with the results in Ref. 16 . 
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FIG. 3. The cross sections in the off-diagonal basis, Eq.(|ll|), at -y/i = 400GeV,u; ma x = lOGeV 
for the e"e + -> tt process: t T tj_ (UD), (DU) and t T t T +t l i l (UU+DD). The suffix "Tree" and 
"SGA" mean the differential cross-section at the tree level and at the one- loop level in the soft 
gluon approximation. It should be noted that DU (UD) component for the e~[e + (e^e + ) process 
is multiplied by 100. 
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FIG. 4. The to max dependence of the cross-sections in the off-diagonal basis at y/s = 400 GeV. 
The DU (UD) component for the e~[e + (e^e + ) process is multiplied by 100. 



t e + 





e t e t 

FIG. 5. The real gluon emission contributions to top quark pair production. 
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FIG. 7. The momentum (unit vectors) configuration of the top and anti-top quarks in the CM 
frame. The momentum of e" (e + ) is in the +z (—z) direction. 
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FIG. 8. The cross sections in the helicity, beamline and off-diagonal bases at \fs = 400 GeV. 
Here we use a "beamline basis", in which the top quark axis is the positron direction in the top 

= 10 GeV. 



rest frame, for both e L e 



and e R e + scattering. For the SGA curves we have used w n 
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Helicity Basis Beamline Basis Off-Diagonal Basis 
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FIG. 9. The cross sections in the helicity, beamline and off-diagonal bases at i/i = 800 GeV. 
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FIG. 10. The cross sections in the helicity, beamline and off-diagonal bases at ^fs = 1500 GeV. 
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